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Recently, the Hamiltonian Control Theory was used in Boreux et at . 20 ldt | to increase the 
dynamic aperture of a ring particle accelerator having a localized thin sextupole magnet. In this 
paper, these results are extended by proving that a simplified version of the obtained general 
control term leads to significant improvements of the dynamic aperture of the uncontrolled 
model. In addition, the dynamics of flat beams based on the same accelerator model can be 
significantly improved by a reduced controlled term applied in only I degree of freedom. 
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1. Introduction 

Hamiltonian Control Theory has been developed in [Vittotl . 120041 : IChandre et all . 120051 ] aiming to improve 
some selected features of a given Hamiltonian system (like to reduce its chaotic behavior), by adding a 
'small control term' which slightly a lters the Ham i ltonia n of the system. 

This technique was adapted in (Boreux et al. . 2010( ] to the case of symplectic maps. In particular, it 
was applied to a four-dimensional (4D) map which models a simplified accelerator ring with sextupole 
nonlinearity, succeeding to increase the stability domain around the nominal circular orbit (the so-called 
dynamic aperture - DA). This procedure allowed the construction of a sequence of control terms of in- 
creasing complexity, that successively improve the dynamics of the original map. Actually, the larger the 
DA became the functional complexity of the control term increased, and the required CPU time for the 



evolu tion of orbits grew. Real accelerators are composed of such simplified maps (see for example [Forest 



1998| |). Thus, we believe that the study of these maps can guide the research for increasing the DA of more 
complicated models. 

The aim of the present paper is to show that also simplified versions of the general control term can 
lead to a significant increase of the DA. For accelerators the addition of a control term can be interpreted 
as the addition of an extra magnet. Since magnetic fields depend only on spatial variables we investigate 
the efficiency of a modified control term obtained by n eglecting the dependence on the momenta of the 
general control term constructed in [Boreux et all l2010f | . 

In many real accelerators the vertical extend of the beam is much smaller than its horizontal size. Such 
cases can be approximated by considering ideal, flat beams of zero height, whose dynamics is described by 
a restriction of the general 4D accelerator model to a 2-dimensional subspace. We also apply our control 
theory to a 2D map model of this type, improving the stability of the flat beam. 



2. The model of a simplified accelerator ring 

As in Bountis fc Tompaidid . 1991 ; Bountis fc Skokos . 2006a ; Boreux et al . 2010 ] we consider a simplified 
accelerator ring with linear frequencies (tunes) q x , q y , having a localized thin sextupole magnet. The 
evolution of a charged particle is modeled by the 4D symplectic map 
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where x\ (x^) denotes the initial deflection from the ideal circular orbit in the horizontal (vertical) direction 
before the particle enters the element, and X2 (x^) is the associated momentum. Primes denote positions 
and momenta after one turn in the ring. The parameters oj\ and 0J2 are related to the accelerator's tunes q x 
and q y by oj\ = 2irq x and UJ2 = 27rq 7/ . In our stu dy, we set q x = 0.61803 and q y = 0.4152 corresponding to a 
regular orbit as shown in [Vrahatis et al\ . \l99lf \. The particle dynamics at the n-th turn, can be described 
by the sequence (xi(n), X2(n), x^{n), £4(n)) n >o, where the (n + l)-th positions and momenta are defined 
as a function of the n-th ones by ([T]). 



3. Theoretical considerations and numerical techniques 

For sake of completeness, let us briefly recall the theoretical framework developed in [Boreux et al\ . \201(t\ ). 
Map (P) naturally decomposes in a integrable part, the rotation by angles co\, L02 in the planes x\,X2 and 
X3, X4 respectively, and a quadratic "perturbation", and can be obtained as the time-1 flow of the following 
Hamiltonian systems 



H(x 1 ,X 2 ,X3,X i ) 
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and V(x\, X2, X3, X4) 



^1 2 
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(2) 



"'^We remark that in the theory we developed q x and q v have been considered as parameters. So, our results can be straight- 
forwardly applied to other cases, remaining valid for tunes satisfying a non-resonant condition (see section below). We also 
note that our method can be easily adapted to the resonant case, but the construction of a new control term is required. 
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Using the notation of the Poisson bracket - which is defined by {H}f := {H, /} = J2j JqJ ~~ Wj^Pj ^ or 
any function f(p, q) - map Q can be written as 

x' = T s {x) = e^e^x, (3) 

where x = (x\, X2, x%, x^) T € M 4 , ( T ) denotes the transpose of a matrix, and the exponential map is defined 

as e^f = E n >o ^f, with {HTf = {IT}*- 1 ({#}/). 

Assuming that the perturbation V is small close to the origin, so that V = o(H), we can construct 
a control map whose generator F is small with respect to V (i.e. it satisfies F = o(V)). Moreover, the 
controlled map 

T ctrl = e W e W e ^ , (4) 

is symplectic and co njugated to a map T*, closer to e^ H ^ than T (for more details the reader is referred to 
' Boreux et qfl . l201Q( |). Assuming furthermore, that a non-resonant condition is satisfied for the particular 



choice of the q x and q y values of map (HJ , the generator F of the control map F } is 



\{V}GV + o{V 2 ) , (5) 



where Q is a "pseudo-inverse" operator, which satisfies Q{1 — e~ n )Q = Q (see Chandre et al . 2005l | for 
more details). 

Even truncating the generator F at order 2 and using 

F2 = \{V}QV, (6) 

as an approximate generator, the control map e^ 2 ^ cannot, in general, be written in a closed analytic form. 
So, we define a truncated control map of order k 

k 



C k (F 2 ) = ^^, (7) 

1=0 

and a truncated controlled map of order k: 

T k {F 2 ) = e ^e^C k (F 2 ) = T s C k (F 2 ) . (8) 

Following (Boreux et al . 2010( |. the Smaller Alignment Index (SALI) method of chaos detection is 
used to determine the regular or chaotic nature of orbits of map ([T|) and its controlled version @. We 
note that an orbit is considered to escape and collide with the accelerator's vacuum chamber if at some 
time n ^i=i a; i( n ) > 10- F° r the evaluation of SALI, the tangent of the studied map is computed and 
employed for following the evolution of two initially linearly independent unit deviation vectors v\ (0) and 
v 2 (0) , defining the index as 

SALI(n) = min{\\vi{n) + £ 2 (n)|[ , \\v\{n) - u 2 (n)||} , (9) 

where II - II denotes the usual Euclidean norm and Vi(n) = rr^r^n , i = 1,2 are vectors of unit norm. 

The behavior of SAL I, and i t s gen e ralization, th e so-called Genera lized Alignment Index (GALI), has 
been studied in detail in [Skokosl . 1200 ll : ISkokos et all l200l 12004! . 12007( 1 . According to these studies, in 4D 
maps the SALI of chaotic orbits tends exponentially to zero as SALI(n) cx e~( ai ~ a2 ^ n (with a±, a 2 being 
the two largest Lyapunov characteristic exponents of the orbit), while it fluctuates around positive values, 
i.e. SALI(n) cx const., for regular orbits. In the case of 2D maps, the SALI tends to zero both for regular 
and chaotic orbits, as SALI(n) oc 1/n 2 and SALI(n) cx e~ ain , respectively. Thus, the completely different 
behaviors of SALI for regular and chaotic orbits allow us to clearly distinguish between the two case both 
for 4D and 2D maps. 

In the following we consider an orbit to be chaotic if SALI(n) < 10 -8 after n iterations. We note that 
in our studies we typically have SALI(n) > 10 -2 for regular orbits. Consequently the chosen threshold 
value (SALI(n) < 10 -8 ) safely guaranties the chaotic nature of orbits. 
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4. The controlled 4D map 

In Boreux et all 120101 ] it was shown that the 4 th order controlled map T^i 7 ^), is a very good choice for 
the controlled system, since it succeeded to considerably increase the DA of the accelerator, keeping also 
the required CPU time for the evolution of large sets of initial conditions, at acceptable levels. 

The computed generator F2 is a complicated function of both the p ositions x-\ , x% and the momenta 



20101 ]. 



X2, Xi of map ([I]), and its specific expression is given in the Appendix of [Boreux et al.\ 

A simple approach f or investigati n g th e global dynamics of 4D accelerator maps, used : 



Bountis fc Skoko; 



2006al : iBoreux et all , hold ], was the computation of the maximal radius of a 4- 
dimensional hypersphere centered at (xi, X2, £3, £4) = (0,0,0,0), containing only regular orbits up to a 
finite but large number of iteration s. This approach pr ovides a reliable indication of the size of the DA. 
Application of this methodology in [Boreux et a/1 12010 ] showed that the controlled map T±{F2) has better 
stability properties than the uncontrolled map ([I]). Nevertheless, one should keep in mind that the center 
of these hyperspheres does not bear any particular physical meaning, as it does not correspond to an 
experimentally realizable beam configuration. 

From a physical point of view it is more meaningful to study the dynamics of the nominal orbit having 
Xi(0) = £3(0) = for different momenta X2, x±. For this purpose, we consider 4-dimensional hyperspheres, 
centered at (0, X2, 0, £4) and compute for different values of X2, £4, the radius R of the largest hypersphere 
which contains only regular orbits. 

The outcome of this investigation is presented in Fig. [TJ where the X2, X4 coordinates of the centers 
of the considered hyperspheres are colored according to the value of radius R for the uncontrolled map 
© (Fig- QI a ))> an d the T±(F2) map ( Fig. CUb)). For both maps the largest value of R is obtained for 
X2 = £4 = (i.e. the value obtained in [Boreux et al. . 2010( ]). while R decreases as the hypershere's center 
is moved away from X2 = X4 = 0. This implies that large values of the momenta introduce instabilities and 
chaotic behavior, which reduce the size of the stability domain. 




Fig. 1. The X2, X4 coordinates of the centers (0, 0, X4) of the 4-dimensional hyperspheres containing only regular orbits 
for (a) the uncontrolled map (Tf]), and (b) the Tn(F2) map. Each point is colored according to the value of the hypersphere 
radius R (with white color corresponding to R = 0), and all orbits are evolved up to 10 4 iterations. 



Nevertheless, the controlled map T4(i ? 2) succeeded to increase the stability region, since R > for larger 
value of X2 and X4, with respect to the original map. In addition, R attains larger values in the central 
region of the (x2, £4) plane for the controlled map T±{F2), being R « 0.6, which suggests a significant 
increase of the stability domain around the nominal orbit of the accelerator. 

As it has already being mentioned, the controlled map T±(F2) has a complicated functional form, 
depending both on spatial variables x\, X3 and on the conjugate momenta X2, X4. From a practical point of 
view, the addition of a control map ([7]) can be thought as the addition of an appropriate magnetic element 
in the ring. On the other hand, the vector potential of a Maxwellian magnetic field which is transverse to 
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the particle motion (i.e. without any longitudinal dependence) is a function of these transverse positions 
and not momenta. It is legitimate then to investigate the efficiency of a modified control term, obtained 
by neglecting the dependenc e of F 2 on momenta . In practice, this means that we set x 2 = X4 = in the 
expression of F 2 obtained in [Boreux et all l20ld |. getting 



F 2 {0 \x 1 ,x 3 ) :=F 2 {x 1 ,0,x 3 ,0) 

1/2 2 
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2 COS (w 2 ) — COS [0J\ + U>2) 

This is a 4 <ft order polynomial in the positions which may be trans formed with s ome y ariable rescaling to 
the vector potential of an octupole magnet with normal symmetry Wiedemann! . 120071 ] 
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where 64 denotes the multi-pole coefficient related to the magnet field strength. Since F 2 does not depend 

on the momenta, the control map 2 , can be explicitly computed. The action of this map on x 2 is given 
by 



e |F 2 } x 2 



x 2 + {F 2 i0) }x 2 + {F2 
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because {i ? 2°^}x 2 does not depend on the momenta and consequently {F^ )> } m x 2 = for all integers 



(0). 



m > 2. Similar relations hold also for the other variables, i.e. {F 2 ^} m Xi 

{FjVyxi = {F 2 
controlled map 



for all integers m > 2 and 



°^} x 3 = 0. Thus, no truncation of the form (J7|) is needed, and consequently the simplified 

(13) 



T, 



(0) 



c 



is symplectic by construction, as a composition of three explicitly known symplectic maps. 

Since map is not explicitly constructed by the Hamiltonian Control Theory, it is intere s ting to 
check its performan ce in controlling map ([T]) and increasing its DA. Following jBountis &: Skokod . l2006al : 
Boreux et al . 2010l ]. the SALI method is used to determine the regular or chaotic nature of orbits of map 



Tq \ In order to directly compare these results with the ones obtained in previous studies, the values of 
log 10 SALI after 10 5 iterations are plotted in Fig. [2], to describe the dynamics in the 2-dimensional subspace 
x 2 (0) = £4(0) = 0, for the uncontrolled map (pQ) (Fig. [2(a)), the T 4 (F 2 ) controlled map © (Fig. 12(b)), and 
the simplified map (]13p (Fig. E(c)). Chaotic orbits are characterized by small SALI values and are 
located in the blue colored domains, regular orbits are colored red, while white regions correspond to orbits 
that escape in less than 10 5 iterations. Fig. [2] shows that both controlled maps, T 4 (F 2 ) andr£ 0) , increase 

the DA of the accelerator. It is remarkable that the simplified controlled map not only increases the 
region of non-escaping orbits, with respect to the original system, but also decreases drastically the number 
of chaotic orbits. 

In order to perform a mo r e global invest i gation of the dynamics of these maps we consider, as was done 



m 



20101 ] . initial conditions inside a 4-dimensional hypersphere 



Bountis Skokosl . I2006al : iBoreux et all 

x 3 = x 4 = 0, and compute the percentages of regular and chaotic orbits as a 



centered at x\ = x 2 



function of the hypersphere radius r (Fig. [3]). From the red curves in Fig. [31 it is observed that the Tq^ 



map significantly increases the domain of regular motion, not only with respect to the original map (black 
curves in Fig. [3j), but also with respect to the controlled map T^{F 2 ) (blue curves in Fig. [3]). In addition, 
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(a) (b) (c) 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



X 1 



Fig. 2. Regions of different SALf values on the (xi, 13) plane of (a) the uncontrolled map {TJ, (b) the 14(^2) controlled map 
|8]>. and (c) the simplified controlled map (|13[) . 16000 uniformly distributed initial conditions in the square (21, £3) G 
[— 1, 1] x [—1, 1], £2(0) = 2:4(0) = are followed for 10 5 iterations, and they are colored according to their final logj q SALI 
value. The white colored regions correspond to orbits that escape in less than 10 5 iterations. 



map has the smallest percentage of chaotic orbits among the studied models, which means that orbits 
of this map either escape or they are regular. 




Fig. 3. The percentages of regular (solid curves) and chaotic (dashed curves) orbits after n = 10 5 iterations of the original 
map U]) (black curves), the 74(^2) controlled map (£8j) (blue curves), and the simplified controlled map (|13|) (red curves), 
in a 4-dimensional hypersphere centered at x\ = X2 = 23 = 24 = 0, as a function of the hypersphere radius r. The largest 
radii at which the percentage of regular orbits is 100%, are also marked. 



Thus, map Tq (|13p controls the original system ([T]) more efficiently than map T4(i ? 2) ©. Additional 
advantages of this map is its simplicity, since it depends only on the spatial variables, and the fact that 
it is symplectic by construction. Another feature, which could be of major practical importance, is that 
the generator (|1U|) of the simplified control map can be considered as the potential induced by a 
magnetic element, since it depends only on spatial variables. 

These results show that the simplified controlled map works better than the controlled map T^{F2) 
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which is also able to enlarge the DA with respect to T$. Then, a natural task would be to compare Tq 
with T ctr i or with T 00 (i ? 2)- We note that we were not able to find a closed analytic form for the latter 
maps due to the cumbersome involved summations. Thus, these comparisons remain open problems which 
we plan to address in the future. However, we remark that the finite number of terms involved in T^i 7 ^) 
possibly allows the attribution of a real magnetic element to this map, while probably this will not be 
feasible for maps T ctr i and T 00 (F2). 



5. Flat beams: the controlled 2D map 

In many particle accelerators the beam is very flat, i.e. its vertical extend is much smaller than the horizontal 
one. A simple, first, approach in investigating the dynamics of a flat beam of the 4D map (pQ) is to neglect 
the vertical coordinate X3 and its corresponding momentum X4, and pass from the initial 4D map to the 
2D map 



X\ 



T 



2D Xl 



(14) 



cos 0J\ — sinwi 

^sinwi coswi J \X2 + x{J \ x ^/ 

whose dynamics was studied for example in [Bountis fc Skokod . l2006bl ]. 

Following for this map, the procedure described in Sect. El truncated control C| D (F 2 2Z3 ) and controlled 
T| D (F 2 2Z) ) maps of order k are constructed, in analogy to Eqs. ([7]) and ([S]), respectively. As before, the 
controlled map is not necessarily symplectic. Since a map is symplectic if its Jacobian matrix A satisfies 

(in its definition domain) the equality A T J A — J = (with J = ( ^ J J being the standard symplec- 
tic constant matrix), the symplectic nature of T% D (F$ D ) can be checked by computing the norm of 
A^JAfc — J, for the Jacobian matrix of the map. The results of this computation are presented in Fig. HI 
for orders k = 4, k = 6, and k = 8, in the region (xi,X2) G [—1, 1] X [—1, 1], and show that T| D (F 2 2Z) ) is 
a good approximation of a symplectic map for k > 6, since < 10~ 6 for a large portion (> 77 %) of 
variable values. As expected, the larger the order k, the closer the map T^ D {F^ ) numerically approaches 
the symplectic condition. 
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Fig. 4. Plot of log 10 D)~, where is the norm of matrix A^JA^ — J, and Af. is the Jacobian of the controlled map T^ D (F2 D ), 
for 16000 uniformly distributed values in the square (xi,X2) £ [— 1, 1] x [—1, 1], for (a) k = 4, (b) k = 6, and (c) k = 8. The 
percentages of orbits with log 10 < —6 are 53%, 77%, and 94% respectively for k — 2, 4, and 6. The color scale corresponds 
to the value of log 10 D^; the smaller the value of log 10 is (blue region), the closer the map is to a symplectic one. 

In order to check whether the controlled map T^ D (F^) with k > 4 increases the DA of the flat beam, 
we use again the SALI to determine the nature of orbits in the (x\,X2) plane, keeping in mind that SALI 
tends to zero both for regular and chaotic orbits of 2D maps, but with time rates which allow the clear 
distinction between the two cases. In particular, many initial conditions are tracked for 10 4 iterations in the 
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(xi, X2) plane, and colored according to their final log 10 SALI value, for the T 2D (Fig. Ela)), the T| D (F 2 2Z) ) 
(Fig. GD(b)), and the Tg D (F2 D ) map. Orbits with SALI < 10 -10 are characterized as chaotic and are colored 
in blue, while the remaining ones are regular. Similarly to Fig. [2] white regions correspond to escaping 
orbits. This distinction is based on the theoretical prediction that for regular orbits SALI oc 1/n 2 = 10 -8 
at n = 10 4 . In Figs. [5^d)-(f) the percentages of regular (blue curves), chaotic (red curves), and escaping 
(green curves) orbits of these three maps are plotted as a function of the radius r of a circle centered at 
the origin x\ = X2 = 0. 




0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 

r r r 



Fig. 5. Upper row: Regions of different SALI values on the (£1,2:2) plane of (a) the uncontrolled 2D map (|14[) , (b) the 
Tg D (F$ D ) controlled map, and (c) the Tg D (F% D ) controlled map. 16000 uniformly distributed initial conditions in the square 
(xi,X2) 6 [—1, 1] x [—1, 1] are followed for 10 4 iterations, and they are colored according to their final log 10 SALI value. The 
white colored regions correspond to orbits that escape in less than 10 iterations. Lower row: The percentages of regular (blue 
curves), chaotic (red curves), and escaping (green curves) orbits after n = 10 4 iterations of the maps of the upper row, in a 
circle centered at x\ = X2 = 0, as a function of its radius r. The largest radii at which the percentage of regular orbits is 100%, 
are marked in panels (d)-(f). Panels in the same columns correspond to the same map. 



In all models the number of chaotic orbits is practically negligible, which means that orbits are either 
regular or escaping. The two controlled maps succeed to increase the DA of the beam, because the domain 
of non-escaping orbits increases (Figs. E)^a)-(c)). Although this domain does not have a cyclical shape, the 
radius of the largest cycle containing only regular orbits increases from r ~ 0.53 for the uncontrolled map 
T 2D CLl), to r « 0.74 for both controlled maps T 6 2D (F 2 2D ) and T 8 2D (F 2 2D ). Since both controlled maps give 
essentially the same results, we conclude that the 6 th order truncation of the controlled map is sufficient 
for controlling the flat beam. 

6. Summary and discussion 



In th is paper, the results of the application of control theory to a simple 4D accelerator map [Boreux et al 



20ld |. are being extended by using simplified versions of the control maps. In the first case, the momentum 



dependence of the control term is artificially removed, so as the control map is a function of the spacial 
coordinates and thus resembles the magnetic vector potential of a multi-pole magnet. In addition, this 
control map is symplectic by construction in contrast to the original one, thus avoiding the associated 
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problems of having to choose an appropriate order of truncation in the Lie representation, for which 
the map satisfies numerically the symplectic condition. The efficiency of the simplified control map is 
remarkable, not only achieving the increase of the DA, but also shows a better performance as compared 
to the complete control map. The remaining important issue regarding the possibility to approximate in 
practice this magnetic field by a multi-pole magnet is still open, and will be addressed in a future study. In 
this respect, and based on the knowledge that control theory can be efficiently used to increase the DA of 
a toy model, the studies can be extended by investigating the possibility of imposing a control map with 
the functional form similar to Eq. (jlip . and compute the associated multi-pole coefficients which achieve 
the best increase in the DA. 

In the second case, the same theory is applied to a 2D version of the map, modeling flat beams, as it is 
the case in electron and positron low emittance rings (i.e. with small beam sizes), where the vertical beam 
size is several orders of magnitude smaller than the horizontal one and thus the dynamics studies can be 
restricted to that plane. This 2D control map presents the same characteristics as the original map, i.e. 
increase of the DA, for a truncation order equal to 6. 

In all our studies, the SALI indicator was mainly used for showing the improveme nt of co n trol i n the 
DA. In future work, we plan to apply the frequency map analysis method (e.g. see Laskarl . Il999l ] and 
references therein) in this map in order to understand the dynamical details of this improvement with 
respect to resonance excitation and diffusion. 
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